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I. 


INTRODUCTION 


The  electron  production  in  the  viscous  layer  over  blunt  bodies  in 
hypersonic  flight  is  dependent  on  the  energy  exchanges  within  the  flow  and 
its  aerodynamic  flight  configuration.  For  a  non-ablating  body  the  reaction 
between  high  temperature  dissociated  air  species  which  requires  the  least 
energy  to  produce  an  electron  is  the  associative  ionization  of  the  nitrogen 
and  oxygen  atoms.  With  increasing  temperatures  collisional  ionization  of  these 
species  begins  to  make  a  contribution  to  the  net  number  of  electrons.  Thus  the 
determination  of  the  electron  density  over  bodies  in  hypersonic  flight  is 
influenced  by  the  vehicle  shape,  the  thermodynamics  of  the  flight  regime,  and 
the  energy  distribution  in  the  various  internal  modes  of  the  gaseous  molecules. 

The  purpose  of  this  report  is  to  predict  the  electron  density  distribu¬ 
tion  over  a  blunt  9°  sphere- cone  body  (Figure  1)  typical  of  the  AFCRL  Trail- 
blazer  II  vehicle.  The  report  is  a  continuation  of  the  study  of  Reference  1. 

In  that  reference  a  discussion  of  the  method  utilized  for  the  prediction  of 
electron  density  was  given  together  with  results  for  several  altitude-velocity 
regimes.  This  report  considers  the  physical  phenomena  (to  be  discussed  sub¬ 
sequently)  of  Importance  to  the  electron  production  and  makes  further  predic¬ 
tions  of  flow  fields  at  different  altitude-velocity  regimes.  The  salient 
features  of  the  predictions  are  given  herein;  the  detail  distributions  of 
velocity,  temperature,  species  concentration  and  in  particular,  the  electron 
density  are  given  in  separate  documents  of  Reference  2. 
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II. 


CHEMICAL  KINETIC  RATE  CONSTANTS 


The  chemical  kinetics  system  for  a  nitrogen- oxygen  atmosphere  is  of 
concern  in  this  study.  The  flight  regime  involves  principally  velocities 
to  17,500  ft/sec.  Since  there  is  no  ablation  of  body  material  the  gas  mixture 
is  composed  of  dissociated  and  ionized  air  species.  For  most  cases  these 

4. 

species  are  O2,  N2,  N,  0,  NO,  NO  and  electrons.  To  cover  tne  conditions 
involving  the  higher  temperature  range  to  15,000°K  and  the  altitude  range  to 

"f*  *f  ^  "f 

100  km  the  additional  atomic  and  molecular  ions  of  Nj  ,  ,  N  ,  and  0  are 

+ 

included.  The  chemical  system  involving  the  neutral  species  and  NO  is  given 

in  reactions  1  to  6  in  Table  I  which  also  contains  additional  reactions  which 

become  important  at  the  higher  temperature  range  of  15,000°K.  This  Table  of 

(3) 

reactions  has  been  suggested  by  Bortner  for  a  nitrogen- oxygen  system 
typical  of  the  shock  layer  of  re-entry  bodies.  At  the  higher  temperatures 
the  associative  ionization  reactions  7  and  8  and  collisional  ionization 
reactions  for  the  neutral  species  is  given  by  9  to  13  become  important. 

Finally,  3fter  the  formation  of  the  positive  ions,  exchange  reactions  redis¬ 
tribute  these  charges  by  means  of  reactions  14  to  20.  The  reaction  rate  con¬ 
stants  have  been  suggested  by  Bortner  in  Reference  (4)  and  (5)  where  he  also 
discusses  their  uncertainties  and  limitations.  Results  in  this  report  are 
obtained  using  reactions  1-6  and  reactions  1-10,  14,  15,  17,  18,  20.  A  minimum 
number  of  reactions  consistent  with  reality  is  utilized. 


In  general,  the  chemical  reactions  are  represented  by  the  equation 


N  N 

£  «ri  Xi  ^  Z  0Ti  Xt 


(1) 


i=l 


t=l 
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where  the  subscript  r  represents  the  particular  reaction,  i  is  the  species 


that  is  being  reacted  and  av^, are  the  stoichiometric  coefficients.  In 
this  report  the  values  of  r  and  i  range  over  the  number  of  reactions  and  the 
number  of  species  respectively.  For  this  representation  the  net  mass  rate  of 


production  of  species  i  per  unit  volume  can  be  immediately  given  by 


R  J  /pc 

=  Mi  rfx  '<Sri)  kfr  B  (  "H7 


-  k  n 

br  j=i 


In  this  equation  is  the  molecular  weight  of  species  i  and  the  subscript  j 
denotes  the  species  and  the  catalytic  bodies,  i.e.,  j  =  1,  2,  ...  N,  N+l,  ...  J. 
The  mass  fraction  of  catalytic  bodies  is  given  in  terms  of  the  N  species  by 


J 

=  M.  L  Z  c,M 

J  i=l  (j-N),i  1  i 


for  j  =  N+l,  ....  J.  The  use  of  equation  (3)  facilitates  the  computation 
problem  for  large  systems  of  chemical  species  and  flow  variables. 

Equation  (2)  is  the  source  term  for  the  mass  generation  of  species  and  is 
a  function  of  the  density  P  of  the  mixture,  the  temperature  T,  and  the  density 
(Pc^)  of  the  pertinent  species.  The  rate  constants  k^r  and  k  are,  of  course, 
functions  of  the  temperature  T. 


III.  VIBRATIONAL  NONEQ JILIBRIUM  EFFECT 


The  excitatior  of  the  molecular  constituents  at  hypersonic  speed  is  a 
complex  phenomenon.  Energy  is  transferred  by  collision  and  is  distributed 
among  the  different  modes  at  different  times  which  are  comparable  to  the 
flight  time,  i.e.,  the  relaxation  times  for  vibrational  and  rotational  modes 
are  different  from  those  for  translation.  Translational  energy  is  rapidly 
absorbed  first  with  the  rotational  mode  requiring  only  a  few  additional  colli¬ 
sions.  In  most  cases  at  the  higher  velocities  and  densities  the  vibrational 
mode  rapidly  becomes  equilibrated  prior  to  the  initiation  of  the  ionization 
reactions.  At  intermediate  values  of  densities  and  temperatures  this  vibra¬ 
tional  equilibrium  time  is  of  the  order  of  the  chemical  time.  In  this  case 
the  vibrational  relaxation  process  has  a  significant  effect  on  the  electron 
production.  Any  change  in  the  electron  density  caused  by  a  lag  in  the  vibra¬ 
tional  mode  depends  on  the  interaction  of  an  increase  in  the  kinetic  tempera¬ 
ture  since  the  heat  capacity  has  not  achieved  its  equiiibium  value  and  a 
decrease  of  the  dissociation  rate  of  the  molecular  species;  thus  the  electrons 
are  decreased  (cr  increased)  since  their  production  depends  mostly  on  the 
dissociation  of  the  molecular  species  except  at  the  very  high  temperatures. 

Ac  the  high  temperature  prevalent  at  the  altitudes  and  flight  velocity  con¬ 
sidered  the  electron  density  is  very  sensitive  to  a  small  change  in  tempera¬ 
ture. 


3.1  Vibrational  Relaxation  Model 

It  has  been  known  for  some  time  that  the  chemical  relaxation  of  the 
flow  after  excitation  by  the  bow  shuck  is  affected  to  an  appreciable  extent  by 
the  state  of  the  vibrational  modes^’^’^.  The  dissociation  of  the  molecular 
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species  depends  largely  on  the  population  of  vibrational  levels  and  the 

probability  of  dissociation  from  each  of  the  individual  levels.  There  is 

a  coupling  between  the  vibrational  relaxation  and  the  dissociation  process 

(7) 

which  has  shown  to  be  important  by  Hammerling,  Teare,  and  Kivel  .  The  effect 
of  dissociation  on  vibrational  relaxation  has  been  studied  by  Treanor  and 

/ON 

Marrone  first'1  '  for  equal  probability  of  dissociation  from  any  vibrational 

level  and  then  for  an  assumed  probability  distribution  which  has  higher  proba- 

(6) 

bility  for  dissociation  from  the  higher  levels  .  By  adjusting  a  parameter 
which  they  callec  the  characteristic  probability  temperature,  good  agreement 
was  found  with  experimental  values  of  dissociation  lag  time  and  rate  constants'''7. 
In  this  model  the  relaxation  of  the  vibrational  modes  occurs  through  a  series  of 
Boltzmann  population  distributions  until  equilibrium  is  reached.  The  principal 
effect  of  this  coupled  vibration-dissociation  is  the  reduction  of  the  dissociation 
rate  constant  from  its  value  at  vibrational  equilibrium.  This,  of  course,  reduces 
the  formation  of  atomic  species.  And  consequently,  the  production  of  electrons 
from  associative  ionization  of  0  and  N  atoms. 


» 


In  the  present  study  the  relaxation  of  the  diatomic  species  09  and  Ng  by 
collision  with  0,  N,  0g  and  Ng  is  examined.  The  average  energy  per  molecule 
obtained  by  summing  over  the  vibrational  levels  of  the  molecules  is  given  by 

*"max 

E  -  Z  (Ei  t)  <4) 

i  =  l 


,  th 


where  n^  and  n  are  the  i  level  and  total  number  densities  respectively. 
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is  the  energy  at  the  i  level.  In  the  consideration  herein  four  collision 
partners  of  0<>,  Nj,  0,  and  N  are  allowed  and  both  molecular  species  0^  and  1^2 
are  vibrationslly  relaxing.  The  relaxation  equation  for  the  energy  as  defined 
by  Equation  (4)  is  then  given  by 


where  n  is  the  number  density  of  the  dissociating  molecules,  Tv  is  the  relaxa¬ 
tion  time,  (dn/dt)f  and  (dn/dt)r  are  the  rate  of  dissociation  and  recombination 
of  the  molecule  and  u  the  velocity  along  the  streamlines.  The  quantities  Ef 
and  Er  are  the  average  vibrational  energy  per  molecule  lost  in  dissociation  or 
gained  in  recombination  respectively.  Ef  is  8  function  of  the  translational 
and  vibrational  temperatures  while  the  energy  lost  to  recombination  is  a 
function  of  the  translational  temperature.  The  energy  at  equilibrium  is 
achieved  when  the  vibrational  temperature  becomes  equal  to  the  translational 
temperature.  Its  difference  from  the  instantaneous  value  of  E  is  the  driving 
function  for  the  relaxing  Landau- Teller  model. 


By  defining  a  quantity  given  by 

JL  -  JL  -  2  I  (6) 

tf  Tv  t  u 

„here  Ty  and  T  ara  the  vibrational  and  translation  temperatures  respectively, 
the  energy  terms  can  be  written  as 
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Er  =  QFU)  f  Ei  exP[“Ei/klJ]  <7b> 

where  the  Q's  are  the  partition  functions  of  the  respective  arguments,  and  as 
before  is  the  energy  of  the  vibrational  levels,  and  U  will  be  given  in  terms 
of  the  ratio  of  the  dissociation  energy  D  and  the  Boltzmann  constant  k. 


The  parameter  U  has  values  that  are  adjustable  and  determines  the  proba¬ 
bility  of  dissociation  for  each  of  the  vibrational  levels.  For  very  large 
values  of  U  there  is  equal  probability  of  dissociation  from  all  vibrational 

levels^.  For  preferential  dissociation  from  the  higher  levels,  values  of 

(3) 

U  =  D/3k  and  D/6k  have  been  used  together  with  a  choice  of  vibrational 

relaxation  times  so  that  the  dissociation  rate  constants  are  in  agreement  with 

_  o 

experiment.  These  values  of  U  have  led  to  a  T  pre- exponential  dependence  of 

the  rate  constant  which  is  too  large  an  exponent  on  the  basis  of  experiment 

data.  The  larger  values  of  U  =  oo  give  a  T  *  dependence  which  is  closer  to  the 

(4) 

value  suggested  by  Bortner  on  the  basis  of  measurement.  However,  the  value 

of  U  =  co  does  not  lead  to  values  of  the  dissociation  lag  time  comparable  to  the 

(9) 

experimental  values  of  Wray  although  the  preferential  model  does  lead  to 

reasonable  values  of  the  dissociation  incubation  time^\  Values  of  the  coupling 

factor  which  is  the  ratio  of  the  actual  dissociation  rate  to  that  at  vibrational 

equilibrium  are  almost  independent  of  temperature  in  the  range  greater  than 

6000°K  for  the  preferential  model  which  has  lower  values  than  the  non-preferential 

(10) 

model.  On  the  basis  of  a  steady  state  analysis  where  the  vibrational  popula¬ 
tion  distribution  which  is  non- Boltzmann,  the  assumption  of  U  -  3T  is  needed  to 
give  a  coupling  factor  with  almost  no  dependence  on  temperature.  Thus  there  is 
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some  choice  in  the  value  of  U  to  be  utilized.  In  this  study,  U  is  kept  as  a 
parameter  with  values  suggested  by  normal  shock  experimental  data  which  can 
be  changed  as  required  by  comparison  with  the  experimental  flight  data. 

In  the  vibrational  model  for  both  the  02  and  N2  molecules  (with  collision 
partners  of  the  principal  air  species  C>2,  N2>  0  and  N) ,  the  relaxation  time  is 
required.  Millikan  and  White have  derived  an  empirical  relaxation  time- 
pressure  relation  as  a  function  of  translational  temperature  for  various  systems 
of  gases.  This  has  been  given  in  the  form 

1  4  /  -i  l\ 

-13  -2  3  /  1  i  \ 

lne  (P Tv)  =  (1.16  x  10  )  fi  $  (  T  -  .015  I  -  18.42  (8) 

where  (p7^)  is  the  product  of  the  gas  pressure  in  atmospheres  and  the  relaxation 
time  Tv  in  sec.,  T  is  °K  in  this  equation,  M  the  reduced  mass  (MiMj)/(Mi  +  M^),# 
is  the  vibrational  constant  h|//k,  and  p  is  in  units  of  atmospheres.  The  range 
of  validity  of  Equation  (8)  extends  o/er  values  of  fi  for  the  air  species 
although  temperature  extrapolation  is  necessary  to  reach  12,000°K  which  prevails 
at  the  shock  front. 

The  molecular  constants  used  for  02  and  N2  are  given  in  Table  II. 

The  vibrational  relaxation  time  is  inversely  proportional  to  the  pressure 
within  the  shock  layer;  this  is  evident  from  Equation  (8).  Values  of  the  time 
are  shown  on  Figure  2  for  the  conditions  from  .003  atm  to  0.1  atm  which  bounds 
the  stagnation  pressure  between  270,000  ft.  down  \.o  200,000  ft.  altitude  for  the 
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flight  velocity  of  17,500  ft/sec.  Characteristic  flow  time  in  the  shock  layer 
is  of  the  order  of  10  “*  sec  at  270,000  ft.  altitude  in  the  stagnation  region. 
This  can  be  compared  with  the  vibrational  relaxation  time  of  Figure  2 .  One 
notes  that  the  characteristic  time  at  the  stagnation  region  can  be  simply 
given  by  the  relation 

tk  ■  [-A  (9) 

char  /  v 

where  y  and  v  are  the  coordinate  and  velocity  normal  to  the  body  surface 
respectively  and  the  integration  limits  are  over  the  fraction  of  the  shock 
layer  where  the  characteristic  time  is  desired. 

The  net  dissociation  and  recombination  from  a  molecular  species  is  the 
sum  of  the  reactions  occurring  at  each  vibrational  level.  In  terms  of  the 
population  distribution  these  rates  are 


i 

max 

n,  i  \ 

kfr 

"  L 

i«=0 

kfr  J  < 

1-max 

I 

-  L 

i=0 

Kr  J 

where  (n^/n)  is  a  Boltzmann  distribution  defined  in  terms  of  the  vibrational 
temperature  T  .  The  summation  is  over  all  the  vibrational  levels.  The  proba- 
bilityv  of  dissociation  from  level  i  is  defined  as 
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I 


P1  ■  Q(Tp)  <“<>  [' VkTF] 


(11) 


i 

and  furthermore,  since  kfr  is  independent  of  Tv<  the  dissociation  rate  constant 
is  given  by 


X  kfr/eq. 


(12) 


The  equation  is  in  terms  of  the  vibrational  equilibrium  rate  constant  and  X  is 
given  by 

Z  exp  [  •  Ej^/kTp] 

Q(TV)  L  exp  [Ei/kuJ  u  ; 

This  function  gives  directly  the  effect  of  the  vibration- dissociation  coupling 
on  the  rate  constant.  In  fact,  this  leads  to  one  effect  on  the  flow  f«.  Id 
electron  density  production  since  the  chemical  rate  constants  are  depressed 
from  the  value  at  vibrational  equilibrium.  A  second  on*',  is  the  increase  of 
kinetic  temperature  at  vibrational  nonequilibrium. 

In  the  above  equation  the  energy  levels  of  the  oscillators  must  be  assigned. 
Experimental  data  can  be  used  to  define  the  energy  levels.  The  partition 
function  Q  can  be  then  simply  written  for  the  harmonic  oscillator  in  the  form 


exp  (-E1QN/kT)  -  1 
exp  (-E10/kT)  -  1 


(14) 


where  N  is  the  number  of  levels  of  the  oscillator,  and  E^q=  hi/  is  the 


vibrational  constant.  The  constants  for  the  0,  and  N2  oscillators  utilized 
in  this  report  are  given  in  Table  II. 

3 . 2  Typical  Results  for  the  Distributions  of  Characteristics  in  t he  Viscous 

Layer 

The  vibrational  nonequilibrium  equations  together  with  the  flow  equations 
and  chemical  source  have  been  coupled  together  and  solved  simultaneously  by  a 
finite  difference  method  utilizing  iteration  to  account  for  the  many  nonlinear- 
ities  of  the  system.  High  altitude  considerations  utilizing  the  shock  slip 
conditions  account  for  the  finite  thickness  of  the  shockweve.  These  equations, 
as  mentioned  before,  have  been  given  in  a  prior  report Consider  the  flow  at 
240,000  ft.  altitude  for  a  body  flight  velocity  of  17,500  ft/sec.  The  Reynolds 
number  at  the  stagnation  region  is  small ^  and  the  entire  flow  region  in  the 
nose  region  is  viscous.  This  can  be  seen  from  the  velocity  curve  on  Figure  3; 
there  is  no  sharp  demarcation  of  slope.  The  velocity  curve  on  Figure  3  also 
shows  the  effect  of  shock  slip.  (The  temperature  was  not  normalized  by  the 
Rankine-Hugoniot  value.)  The  species  distributions  are  denoted  by  t.he  corres¬ 
ponding  name  on  Figures  3  and  4.  Peak  values  of  both  oxygen  and  nitrogen  atoms 
occur  at  about  5  mm  distance  away  from  the  body  whereas  nitric  oxide  peaks  at 
about  2  mm.  The  peak  temperature  occurs  near  the  edge  of  the  layer  where  the 
dissociation  of  molecular  oxygen  is  also  maximum.  Thus  convection  of  the 
dissociated  species  from  the  edge  together  with  diffusion  leads  to  the  distri¬ 
bution  shown  on  Figure  4.  Correspondingly,  the  electron  production  maximizes 
near  the  outer  edge  of  the  layer  whereas  the  electron  density  peaks  nearer  the 
body  surface  (at  about  2.6  mm  from  the  surface  of  Figure  5).  There  are  three 
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temperatures  in  the  viscous  layer  and  these  are  the  kinetic  temperature  and 
the  vibrational  temperatures  of  the  oxygen  and  nitrogen  molecules.  The 
distributions  of  these  temperatures  are  shown  on  Figure  6.  The  vibrational 
temperature  of  the  oxygen  molecules  rises  rapidly  near  the  shock  front  and 
approaches  the  kinetic  temperature  near  the  wall.  Of  particular  significance 
here  is  the  fact  that  vibrational  nonequilibrium  occurs  through  most  of  the 
shock  layer  thickness  so  that  the  reaction  rate  constants  do  not  have  their 
full  equilibrium  values  and  there  is  r  correspondingly  higher  kinetic  tempera¬ 
ture  than  for  the  vibrational  equilibrium  case.  The  electron  density  production 
is  changed  in  opposite  directions  by  these  two  effects;  the  effect  that  domi¬ 
nates  the  electron  density  ic  dependent  on  the  altitude.  The  change  in  electron 
density  with  altitude  is  shown  on  Figure  7.  It  is  important  to  note  that  the 
vibrational  nonequilibrium  effects  on  electron  density  is  contained  principally 

in  the  altitudes  shown  therein.  This  is  an  example  of  typical  results; 

(2) 

additional  ones  are  contained  in  the  report  referenced  previously  and  in  the 
1970  report^  .  Figure  7  shows  that  between  the  altitudes  shown  the  electron 
density  for  the  vibrationally  nonequilibrium  flow  is  higher  than  that  for 
vibrational  equilibrium  with  a  cross-over  at  about  270,000  ft.  altitude.  The 
principal  phenomena  controlling  these  distributions  are  the  dissociation  of  ths 
oxygen  molecule  and  the  magnitude  of  the  kinetic  temperature  which  is  dependent 
on  the  degree  of  vibrational  equilibrium.  Near  the  altitude  of  270,000  ft. 
the  dissociation  rate  of  oxygen  molecules  for  vibrational  equilibrium  is 
higher  than  that  for  vibrationally  nonequilibrium  despite  the  higher  kinetic 
temperature. 
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IV.  EFFECT  OF  VARIOUS  PHENOMENA  ON  THE  ELECTRON  DENSITY  PRODUCTION 


There  are  a  number  of  factors  that  can  influence  the  magnitude  of  the 
electron  density.  Considerations  are  given  to  them  in  this  Section. 

4. 1  Wall  Catalyticity 

The  wall  temperature  at  the  high  altitudes  considered  in  this  Report  is 
low  and  it  is  expected  that  the  dissociated  species  recombine  in  the  gas  phrse 
to  their  molecular  constituents.  Experimental  measurements  of  gas  radiation 
appear  to  indicate  that  nitric  oxide  has  a  non- zero  value  near  the  surface  to 
account  for  the  magnitude  observed.  If  this  is  the  case  then  the  electron 
density  is  increased  from  that  of  a  fully  catalytic  surface.  The  non- catalytic 
condition  for  species  is  given  by  the  vanishing  of  the  mass  flux 

i  =  0  Cl  5) 

where  i  denotes  the  species.  Since  there  is  no  mass  transfer  at  the  body  sur¬ 
face  Equation  (15)  implies  a  similar  condition  on  the  diffusion  flux. 

Consider  the  non- catalytic  wall  condition  for  the  nitric  oxide  con¬ 
centration  only.  There  are  no  atomic  species  at  the  wall.  The  electron  density 
distribution  has  been  obtained  with  this  condition  for  the  altitude- velocity  of 
270,000  ft.  -  17,500  ft/sec  for  vibrational  equilibrium.  A  9.17.  increase  in 
peak  electron  density  over  the  case  with  a  fully  catalytic  wall  condition  was 
obtained.  A  comparison  of  electron  density  changes  due  to  this  condition 
together  with  other  effects  are  shown  on  Table  III. 
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It  would  be  of  interest  to  compare  the  electron  density  for  the  two 
extremes  of  wall  conditions  and  these  are  a  fully  catalytic  one  and  a  non- 
catalytic  one  for  all  species.  Of  course,  at  the  low  wall  temperature  con¬ 
cerned,  it  is  expected  that  the  former  condition  would  be  more  typical  for 
heat  sink  body  material.  For  the  non- catalytic  wall  for  all  the  species  there 
are  non- zero  values  of  the  atonic  species  in  addition  to  those  of  nitric  oxide. 
The  results  for  the  240,000  ft.  altitude  and  17,500  ft/sec  velocity  case 
(vibrational  equilibrium)  show  that  the  nitric  oxide  and  its  ion,  and  the 
nitrogen  atom  concentrations  (mass  fraction)  still  peak  within  the  layer  but 
the  other  species  do  not.  However,  since  the  mass  density  increases  toward 
the  wall  the  electron  density  achieves  its  maximum  there.  The  electron  density 
distributions  for  both  wall  conditions  are  compared  on  Figure  8.  The  peak 
values  differ  by  an  order  of  magnitude  as  expected.  As  pointed  out  before  it 
is  likely  that  the  wall  condition  is  non- catalytic  for  all  species.  Air  Force 
Cambridge  Research  Labs,  experimental  flight  data  jjRotman^^,  Hayes 
appear  to  indicate  that  the  electron  density  increases  away  from  the  surface  at 
least  for  streamwise  points  along  the  body  away  from  the  stagnation  region.  It 
is  noted  that  the  fully  catalytic  wall  conditions  for  atomic  species  also  lead 
to  the  same  type  of  variation.  Thus  the  fully  non-catalytic  wall  condition 
would  certainly  give  the  maximum  electron  density  obtainable  as  far  as  recom¬ 
bination  wall  effects  are  concerned. 

4.2  Slip  Effects  at  the  Wall 

Gas  molecules  exchange  energy  and  momentum  with  the  surface  of  the  body 
by  means  of  collision.  Since  there  are  velocity  and  temperature  differences 
between  the  gas  and  the  body  surface,  momentum  and  energy  are  transferred  bet- 
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ween  them.  When  the  gas  molecules  do  not  attain  the  momentum  and  energy  corres¬ 
ponding  to  the  conditions  at  the  boundary  in  a  single  collison  then  there 
exists  a  velocity  slip  and  a  lack  of  complete  accommodation  between  the  gas 
and  the  wall.  Thus  velocity  and  temperature  jumps  occur  at  the  solid  boundary 
over  a  distance  of  the  mean  free  path  of  the  gas.  In  a  gas  with  high  density 
(low  altitude  flight)  these  effects  are  not  important,  they  only  became 
important  in  the  high  altitude  regime.  The  effects  of  slip  and  the  lack  of 
accommodation  result  in  discontinuities  of  the  tangential  velocity  and  the 


temperature  and  are  given  for  the  velocity  u  and  temperature  T  by  relations 


(14,15) 


u(0) 


^  I 


(ft). 


(16a) 


T(0) 


2-Cf 


w 


2  y  X_ 

y+ 1  pr 


d  t 
d  y 


(16b) 


where  the  coefficients  multiplying  the  derivatives  on  the  right  side  of  Equations 
(16a,  16b)  are  known  as  the  slip  coefficient  and  the  temperature  jump  coefficient 
respectively.  In  addition,  L  denote  the  local  mean  free  path  of  the  gas  and 
becomes  small  for  high  density,  f  the  isentropic  constant,  O  is  the  fraction  of 
diffusely  reflected  molecules  whose  average  tangential  velocity  is  zero,  &  is 
the  accommodation  coefficient  (temperature),  and  Pr  is  the  local  Prandtl  number. 
It  is  seen  that  the  jumps  are  proportional  to  the  mean  free  path  of  the  gas. 
Equations  (16a  and  b)  are  derived  from  serai-macroscopic  agreement  usually  for  a 
monatomic  gas;  Equation  (16b)  is  written  in  the  form  valid  for  a  polyatomic  gas. 
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Recent  examination^^  of  the  temperature  jump  in  a  polyatomic  gas  shows  that 

the  values  of  the  temperature  jump  are  always  slightly  less  than  those  for 

monatomic  gases.  The  two  coefficients  O'  and  Ot  are  prescribed  constants  and 

(14) 

are  close  to  unity.  For  example,  Oi  vary  from  .87  to  .97  for  aluminum. 

The  values  were  taken  to  be  unity  for  the  numerical  results  given  in  this 
report.  It  is  seen  from  Equations  (16a  and  b)  that  both  the  velocity  deriva¬ 
tives  are  positive  for  a  cool  wall.  This  fact  implies  that  there  are  positive 
jumps  in  both  velocity  and  temperature  which  are  inversely  proportional  to  the 
ambient  air  density. 

An  evaluation  of  this  effect  has  been  made  for  the  (270,000  ft;17,500  ft' 
sec)  condition.  At  the  wall  the  velocity  ratio  has  jumped  to  .03147  from  the 
zero  value.  The  temperature  in  the  gas  at  the  wall  has  increased  by  a  factor 
of  1.9  of  the  body  temperature.  This  has  the  consequence  of  increasing  the 
magnitude  of  the  temperature  distribution  (Figure  9)  and  therefore  the  con¬ 
centration  of  atomic  species.  The  net  increase  in  the  peak  number  density  of 
electrons  for  this  case  is  only  1.77,  and  is  3mall.  Although  the  temperature  near 
the  body  surface  has  increased  sufficiently  to  noticeably  influence  the  ioniza¬ 
tion  rate  there,  most  of  the  electrons  are  produced  principally  near  the  outer 
two- thirds  part  of  the  layer  and  hence  the  peak  electrons  are  not  noticeably 
changed.  Near  the  wall  the  concentration  of  nitric  oxide  ion  is  increased  by 
the  higher  temperature  due  to  the  jump;  however,  the  electron  number  density  is 
decreased  by  2L7.  due  to  the  fact  that  the  density  is  decreased  by  the  temperature 
jump. 
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Since  the  jumps  of  velocity  and  temperature  are  proportional  to  their 
derivatives  and  mean  free  path  their  variation  around  the  body  requires 
further  evaluation. 

4 . 3  Effect  of  Diffusion  Coefficient 

Diffusion  has  a  considerable  influence  in  distributing  the  electrons 

produced  in  the  high  temperature  region  near  the  shock  front.  This  fact  has 

been  pointed  out  in  a  previous  paragraph.  The  results  given  in  this  report 

have  been  obtained  for  a  constant  Lewis  number  (=  PD^jCp/k  where  p  is  the  mass 

density,  the  binary  diffusion  coefficient,  Cp  the  mixture  specific  heat  at 

constant  pressure,  and  k  the  coefficient  of  heat  conductivity)  in  which  P, 

C  ,  and  k  are  computed  for  each  point  as  a  function  of  the  mixture  concentration 
P 

and  temperature.  Comparisons  have  shown  that  the  assumption  of  a  constant  Lewis 
number  is  a  valid  representation  for  a  multi- component  diffusion  providing  there 
are  no  large  disparate  species  masses.  Thus  dissociated  air  species  satisfy 
this  requirement;  an  example  in  which  it  is  not  satisfied  is  a  hydrogen- air 
mixture . 

A  Lewis  number  of  1.4  has  been  utilized  in  all  the  results  considered  in 
this  report  except  as  noted  differently.  To  demonstrate  the  effect  of  diffusion 
numerical  results  for  the  case  (270,000  ft; 17, 500  ft/sec)  have  been  obtained 
with  the  Lewis  number  equaled  to  1.6.  The  peak  electron  density  decreases  by 
17.57,  due  to  increasing  the  diffusion  coefficient.  That  is,  there  is  a  higher 
spread  of  the  electrons  produced. 


Effect  of  Transverse  Steo  Size  on  Electron  Density  Near  Wall 


Since  there  are  flight  measurements  which  are  close  to  the  surface  of 
the  body  the  nonlinear  variation  in  that  region  is  of  particular  interest. 
This  point  was  investigated  by  decreasing  the  uniform  step  size  and  by  the 
use  of  a  coordinate  transformation  which  created  a  nonuniform  step  size 
normal  to  the  wall  with  a  smaller  step  size  near  the  wall  as  required.  With 
use  of  the  transformation  the  second  order  accuracy  in  the  finite  difference 
scheme  was  retained.  The  independent  variable  Tf  transverse  to  the  wall  is 
transformed  to  the  £  coordinate  by  the  function 


t  =  Hr 


tanh  L/V 


tanh  Ct 


where  a  is  a  constant  to  be  assigned.  This  constant  is  used  to  adjust  the  step 
size  near  the  wall.  Smaller  values  of  a  will  lead  to  higher  concentration  of 
increments  near  the  wall.  The  equations  of  motion  as  given  in  Reference  (1) 
were  transformed  to  the  new  coordinate  and  modifications  were  made  accordingly. 

An  example  of  the  electron  density  variation  for  different  step  sizes  is 
shown  on  Figure  10. 

4 . 5  Effect  of  Chemical  Rate  Constants  on  the  Formation  of  Electrons 

It  was  pointed  out  in  Section  II  that  at  the  higher  temperatures  and  low 
density  additional  ionization  mechanisms  may  assume  some  importance  in  the 
production  of  electrons.  Consideration  is  given  in  this  section  to  the  behavior 
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of  the  electron  density  for  different  chemical  rate  constants.  A  first 
example  is  the  flight  condition  at  the  highest  altitude  (270,000  ft)  which 
has  a  translational  temperature  near  the  shock  front  of  about  9630°K.  It  is 
noted  the  all  shock  layer  temperatures  for  flight  conditions  between  the 
altitude  200,000  ft.  and  270,000  ft.  for  the  velocity  17,500  ft/sec  are 
below  11,500°K.  Most  values  of  the  temperatures  are  below  10,000°K  except 
for  the  one  case  at  200,000  ft.  altitude.  Results  for  the  (270,000  ftjl.7,500 
ft/sec)  case  with  the  additional  ionization  mechanisms  given  by  re<_  tions  1-10, 
14,  15,  17,  18,  and  20  of  Table  I  were  obtained.  With  these  reactions  the 

*4*  -f  *f 

species  of  Ng  ,  0^  ,  0  ,  N  are  added  to  the  other  air  species.  The  peak 

electron  density  is  increased  by  a  factor  of  5  with  the  additional  reactions 

and  species.  A  comparison  of  the  distribution  of  electrons  between  the  5 

species  and  10  species  cases  is  given  on  Figure  11.  There  are  now  five  positive 

+  4* 

ions  generated  within  the  flow.  Of  particular  interest  is  that  NO  and  02  are 
almost  equal  and  are  larger  by  one  to  two  orders  of  magnitude  than  the  N+  and  0+ 
ions.  This  difference  is  shown  on  Figure  12.  Although  the  electron  density  has 
not  changed  appreciably  the  dominant  ion  now  is  0^  .  Most  of  the  concentration 
of  (>2^  ion  is  produced  near  the  shock  front  and  the  production  rate  at  this 
point  is  higher  than  that  of  N0+  by  an  order  of  magnitude. 

A  study  was  also  made  on  the  effect  of  changing  the  various  reaction  rate 
constants  on  the  production  of  electrons  for  6  species.  The  dissociation  rate 
of  oxygen  molecules  (reaction  1)  was  increased  by  an  order  of  magnitude  for  the 
(240,000  ft.  altj  17,500  ft/sec)  vibrational  equilibrium  case  and  resulted  in 
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producing  more  atoms  as  expected.  The  peak  temperature  was  decreased  by 
about  12%  and  resulted  in  a  34%  decrease  in  the  peak  electron  density.  It 
may  have  been  expected  that  the  increase  in  oxygen  atoms  would  result  in  a 
higher  electron  density  but  as  it  is  seen  the  dominant  factor  here  is  the 
temperature.  Thus  the  flow  field  variables  clearly  can  not  be  uncoupled  from 
the  species  production. 

In  order  to  determine  the  effect  of  the  deionization  rate  of  reaction  6 
its  rate  constant  was  lowered  by  one  order  from  the  value  given  in  Table  I. 

This  resulted  in  a  small  increase  in  the  electron  density  of  2.5%  for  the 
vibrational  equilibrium  (240,000  ft;17,500  ft/sec)  case.  Thus  the  electron 
density  distribution  has  not  been  changed  appreciably  by  this  rate  constant. 

The  electron  density  distributions  are  shown  on  Figure  11. 

4.6  Variation  of  Electron  Density  Along  the  Body 

The  entire  distribution  of  electron  density  for  vibrational  equilibrium 
(270,000  ftjl7,500  ft  sec)  over  the  body  has  been  considered  utilizing  the 
method  reported  in  Reference  (1)  for  six  species  and  the  chemical  reactions  1-6 
on  Table  I.  Two  sets  of  typical  variations  of  electron  density,  velocity  ratio, 
and  temperature  ratio  are  shown  on  Figure  (13)  and  (14)  for  the  two  locations 
at  6.6  cm  and  76.45  cm  along  the  surface  measured  from  the  stagnation  point. 

It  is  seen  that  the  temperature  peak  develops  within  the  viscous  layer  back 
along  the  body  as  expected  for  blunt  conical  bodies.  The  peak  electron  density 
along  the  body  is  shown  on  Figure  (15);  the  maximum  occurs  near  the  stagnation 
region  and  decays  toward  the  local  equilibrium  value. 
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V. 


CONCLUDING  REMARKS 


The  electron  density  for  the  viscous  flow  over  a  body  in  !  ypersonic 
flight  and  high  altitudes  has  been  studied  in  the  report.  A  number  of  phenomena 
which  can  influence  the  production  of  electrons  has  been  quantitatively  evalu¬ 
ated.  These  include  the  effect  of  vibrational  relaxation  of  the  oxygen  and 
nitrogen  molecules,  wall  catalyticity ,  slip  effects  at  the  wall,  magnitude  of 
the  diffusion  coefficients,  magnitude  of  the  chemical  kinetic  rate  constants, 
chemical  kinetic  mechanisms  at  high  temperature,  and  transverse  step  size  within 
the  flow  field.  The  principal  conclusions  can  be  summarized  as  follows: 

1.  Vibrational  nonequilibrium  effects  cause  higher  electron  density 
between  200,000  ft.  and  slightly  below  270,000  ft.  altitudes  for 
the  flight  velocity  of  17,500  ft /sec. 

2.  The  non- catalytic  wall  condition  for  nitric  oxide  increases  the 
electron  density.  Non- catalytic  condition  for  all  species  causes 
the  largest  increase  of  electron  density  over  that  of  a  ful., 
catalytic  wall. 

3.  Slip  effects  at  the  wall  due  to  the  fact  that  the  energy  and 
momentum  transferred  between  the  gas  and  wall  are  not  accomplished 
in  a  single  collision  are  localized  near  the  wall.  There  is  an 
increase  in  the  N0+  concentration  but  a  large  decrease  in  electron 
density  near  the  wall  due  to  the  temperature  jump.  Further  study 

is  required  to  determine  its  effect  along  the  body  and  its  effect  on 
the  additional  positive  ions  at  the  higher  temperature. 
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4.  Larger  diffusion  coefficient  lowers  the  peak  electron  density  and 
distributes  the  electrons  to  all  parts  of  the  layer.  V'ithin  the 
spread  of  Lewis  number  of  1.4  and  1.6  there  is  a  17.5%  decrease. 

5.  The  nonlinear  variation  of  electron  density  at  the  wall  has  been 
obtained  for  small  step  sizes. 

6.  Chemical  kinetic  mechanisms  and  rates  for  the  inclusion  of  the 
nitrogen  and  oxygen  molecular  and  atomic  ions  which  are  present  at 
the  higher  temperatures  have  been  considered  for  a  total  of  10  species 
and  15  reactions.  One  case  of  (270,000  ft.  alt;17,500  ft/sec)  at  the 
stagnation  region  has  been  considered  and  the  peak  electron  density 
has  been  increased  by  a  factor  of  5  from  that  of  the  6  species  model. 

7.  The  effect  of  increasing  the  dissociation  rate  of  oxygen  decreases  the 
peak  electron  density  even  though  the  atomic  oxygen  has  been  increased. 

The  reason  is  that  the  temperature  has  been  lowered. 

8.  An  order  of  magnitude  change  of  the  deionization  rate  of  the  nitric 
oxide  ion  has  only  a  small  effect  on  the  electron  density. 

These  results  have  been  obtained  from  the  study  of  a  number  cf  cases 
limited  to  the  stagnation  region.  An  evaluation  of  the  electron  density  at 
other  conditions  at  the  stagnation  region  and  at  points  around  the  body  will 
add  further  information. 

Detailed  numerical  results  are  given  in  separate  documents  of  Reference  (2). 
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CHEMICAL  SYSTEM  OF  AIR  REACTIONS  AND  RATES  (Continued) 


L  M 


TABLE  I 

CHEMICAL  SYSTEM  OF  AIR  REACTIONS  AND  RATES  (Continued) 


Third  Body  Efficiencies 


Species 


'o  25  1  1.1  1 

N  25  3.24  1.1  1 

2 

0  105  1  50  1 

N  I  21.6  50  1 

NO  1  1  50  1 


f 
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TABLE  II 


MOLECULAR  CONSTANTS 


Species 

E1q  =  hV  ERGS 

Dissociation 

Energy,  D/k°K 

Levels* 

(N) 

°2 

-13 

3.1409655  x  10 

59400 

26 

N2 

4.5897248  x  10"13 

113200 

34 

J.u  O 

k  =  1.38  x  10  Ergs/  K,  Boltzmann  Constant 

-27 

h  =  6.625  x  10  Ergs  sec,  Planck's  Constant 

V  =  frequency 

*  Cut-off  Harmonic  Oscillator 
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TABLE  III 

PEAK  ELECTRON  DENSITY  CHANGE  FOR  VARIOUS  CONDITIONS 

Condition  Percent  Change  in  Peak 

Electron  Density 

Nitric  Oxide  -  Noncatalytic  +9.1 


Non- catalytic  -  All  Species  +600 


Wall  Slip  +1.7  (Electron  density  near  the 

wall  is  decreased  by  2\°L) 


Additional  Ionized  Species  +500 


Increase  Oxygen  Dissociation  -34. 

Rate  -  One  Order 


Decrease  Deionization  Rate  +2.5 

(Reaction  6)  by  One  Order 
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Figure  3.  Vibrational  Nonequilibrium  Distributions  of  Velocity, 

Temperature,  and  Species  Concentration  in  Mass  Fraction 
at  Stagnation  Region,  2^0,000  ft;  17,500  ft/sec 
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Figure  5. 


Electron  Density  (CC  *)  for  Vibrational  Nonequilibrium  at  the 
Stagnation  Region,  240,000;  17,500  ft/sec 
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Figure  7.  Peak  Electron  Density  Variation  with  Altitude  at  the 
Stagnation  Region  for  the  Velocity  of  17,500  ft /sec 
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Figure  8.  Electron  Density  Profiles  for  Catalytic  and  Non-Catalytic  Walls 
at  the  Stagnation  Region,  240,000  ft;  17,500  ft /sec 
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Figure  10.  Effect  of  Step  Size  on  the  Distribution  of  Electron  Density  at 

the  Stagnation  Region  for  Vibrational  Nonequilibrium,  240,000  ft 
17,500  ft 'sec 
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Electron  Density  Distributions  at  the  Stagnation  Region  for 
6  Species  (N0+)  and  10  Species  (N0+,  Ng  ,  0^  >  N  ,  0+)  at 
270,000  ft  and  17,500  ft/sec 
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Figure  11. 
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Figure  13.  Temperature,  Velocity,  and  Electron  Density  Profiles  for  the 
Sphere  Cone  Trailblazer  Body  at  x  =  .066  meter,  270,000  ft; 
17,500  ft/sec 
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Figure  14.  Temperature,  Velocity,  and  Electron  Density  Profiles  for  the 
Sphere  Cone  Trailblazer  Body  at  x  =  .066  meter,  270,000  ft; 
17,500  ft/sec 
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